/* Integration points (quadratures)
   Jan Stebel
*/

#ifndef _QUADDISCPOINT_H__
#define _QUADDISCPOINT_H__


#include "quadpoint.h"


/** Quadratic discontinuous interpolation */
template<unsigned int dim>
class QuadDiscPoint : public QuadPoint<dim>
{
protected:
	using IntPoint<dim>::imesh;
	bool isLinear() { return false; };
	bool isDiscontinuous() { return true; };

public:
	QuadDiscPoint() {};
	QuadDiscPoint(Mesh *mesh, const Quadrature<dim> *q) { QuadPoint<dim>::init(mesh, q); };
	int mapDof(int e, int n);
	int mapDof(int e, int n, int p);
	inline int getMeshNdof() { return imesh->getNvolumes()*10 + imesh->getNfaces()*6 + imesh->getNedges()*3; }
};










#endif
